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1.  INTRODUCTION  AND  PROJECT  OBJECTIVES 


Nonlinear  optics,  the  study  of  matter  interactions  with  intense  electromagnetic 
fields,  has  uncovered  surprising  and  useful  physical  effects,  such  as  stimulated 
scattering  of  light  by  which  an  intense  monolithic  beam  can  be  converted  to  intense 
coherent  beams  of  different  wavelengths.  Recently,  interesting  processes  involving 
nonlinear  optical  image-bearing  beams  have  been  discovered,  such  as  optical-beam 
phase  conjugation.  Optical-beam  phase  conjugation  is  the  name  given  to  any  process 
which  generates,  in  real  time,  the  time-reversed  replica  of  a  complex,  image-bearing, 
optical  beam,  or  which  generates  other  related  reflected  beams,  both  monochromatic 
and  polychromatic.  Phase  conjugation  by  optical  four-wave  mixing  was  conceived 
earlier  in  this  project,  as  were  such  applications  of  phase  conjugation  as  the  corrections 
image  aberrations,  brightness  enhancement  of  laser  outputs,  edge-enhancements  of 
images,  automatic  steering  of  beams,  and  optical  spectroscopy  by  Raman-induced 
phase  conjugation  (RIPC).  Such  nonlinear  optical  techniques  are  currently  being 
extended  to  perform  many  forms  of  beam  processing,  sorting,  and  routing  on 
picosecond  time  scales.  This  project  aims  at  exploring  and  developing  these  and  other 
new  wave-mixing  processes,  and  the  necessary  nonlinear  materials.  Other  novel 
electromagnetic  scattering  processes  which  occur  at  intense  optical  fields  available 
from  lasers  are  also  explored  for  their  scientific  and  device  interest.  The  approaches 
taken  in  this  project  are  both  experimental  and  theoretical. 


2.  MAJOR  ACCOMPLISHMENTS:  12/1/93  to  4/30/97 

The  major  accomplishments  of  this  grant  from  the  last  date  covered  by  the  previous  Tech¬ 
nical  Report  (8/3 1/96)  to  the  final  date  of  the  grant  ( 4/30/97)  are  described  briefly  below  and 
more  completely  in  the  numbered  references  and  the  Appendices  I  and  n  attached  to  this  report. 
The  major  accomplishments  from  the  beginning  of  this  grant  up  to  8/31/96  are  copied  from  the 
three  yearly  Technical  Reports  proceeding.  These  give  numbered  references  to  complete  descrip¬ 
tions  of  the  work,  which  references  appear  in  Sections  3  and  4  of  this  report. 

2.1  A  new  interferometric  technique  to  measure  amplitude  and 
phase  of  nonlinear  tensor  elements. 

We  have  demonstrated  an  all-optical  method  that  allows  the  measurement  of  the  amplitude 
and  the  phase  of  all  independent  elements  of  the  degenerate  four- wave-mixing  tensor  of  an  isotro¬ 
pic  thin  film  on  a  thick  substrate.  The  method  employs  the  interference  of  the  degenerate  four- 
wave-mixing  signal  from  the  substrate  with  that  from  the  film,  normalized  to  other  signals  gener¬ 
ated  by  the  same  incident  beams.  Our  observation  of  a  negative  real  part  in  every  element  of  % 
at  768  nm  of  supports  the  existence  of  a  two-photon  resonance,  suggested  by  others  to  exist 

near  1.3  |im .  This  research  is  described  in  detail  in  Ref.  3.21  which  is  attached  as  Appendix  I. 

2.2  Observation  of  the  mobility  of  a  large  polaron 

We  completed  a  series  of  researches  on  the  nature  of  photoconduction  in  n-type  cubic 
Bij2SiO20-  As  a  result  it  became  clear  that  the  conduction  electrons  in  this  crystal  constitute  the 
strongest  physical  example  of  a  long-postulated  theoretical  entity,  a  “large  polaron.”  The  argu¬ 
ments  leading  to  this  conclusion  are  given  in  Pub.  3.22  which  is  attached  as  Appendix  II.  Our 
results  were  also  presented  at  the  conferences  4.12  and  4.13  below. 

2.3  The  second  electro-optic  coefficient  of  benzene 

We  completed  the  first  measurement  of  the  “second”  electro-optic  coeffecient  of  liquid 
benzene  which  is  often  used  as  a  nonlinear  standard.  The  results  are  described  in  Pub.  3.20  and  the 
conference  talk  4. 14. 


MAJOR  ACCOMPLISHMENTS:  12/1/93  TO  10/31/94 


The  major  accomplishments  of  this  grant  period  were  as  follows. 

Atoms  for  logic:  an  all-optical  image  correlator. 

During  the  current  period  we  achieved  one  long-standing  goal:  to  demonstrate  an  optical 
correlator  that  is  fest  (can  correlate  more  than  one  million  images  per  second),  high  capacity 
(more  than  50,000  pixels  per  square  centimeter  of  nonlinear  (vapor)  material),  and  requires  low 
power  (we  need  only  3  mW  for  1  microsecond  to  correlate  a  36  pixel  pattern  with  a  pattern  of 
280  pixels).  The  nonlinear  medium  was  a  simple  glass  cell  containing  cesium  vapor  at  57°C.  The 
images  were  digital  arrays  which  allow  quantitative  assessment.  We  believe  that  this  correlator 
has/the  best  overall  figure-of-merit  of  any  all-optical  image  correlator  yet  demonstrated,  and 
stands  as  the  standard  for  all-optical  correlators. 

Details  of  our  device  and  its  performance  are  given  in  papers  3.4  and  3.5  listed  below. 
Developments  in  preceding  technology  of  nonlinear  vapor  cells  are  also  given  in  papers  j.8  and 
3.11  below. 


Nonlinear  optical  susceptibilities  of  polymer  films  and  solutions. 

We  continued  our  collaborative  effort  to  develop  better  nonlinear  optical  polymers  with 
L.R.  Dalton  and  colleagues  of  the  USC  chemistry  department,  and  with  C.W.  Spangler  and 
coUeagues  at  Northern  Illinois  University.  We  concentrated  our  efforts  on  recently-developed 
optical  polymers  that  have  greatly  improved  thermal  stability  and  solubility,  especially  bis-thienyl 
polyenes.  Our  results  are  summarized  in  papers  3.1,  3.3  and  3.7  listed  below. 


Nonlinear  optics  of  fullerene  films  and  solutions. 

We  continued  our  investigations  of  the  physics  of  nonlinearities  in  films  and  solutions  of 
the  cage-molecules  C50  ^70  called  "fullerenes".  We  made  an  improved  experimental 

determination  of  the  large  excited  state  absorption  cross-section  of  the  Cgo  molecule.  This  cross- 
section  is  crucial  in  the  design  of  optical  limiters.  Details  are  given  in  paper  3.2  below. 

We  made  the  first  measurement  of  the  dependence  on  wavelength  fi'om  740  to  880  run,  of 
the  nonlinear  susceptibility  tensor  elements  of  the  C50  molecule.  Both  magmtude  and  phase  of 
these  elements  were  determined  as  we  described  in  paper  3.5  below. 

We  also  performed  a  simple  model  calculation  of  the  hyperpolarizability  of  the  Cgo 
molecule.  The  results  are  summarized  in  paper  3.10  below. 

'  New  techniques  for  measurement  of  nonlinear  optical  coefficients. 

We  developed  a  simple  method  by  which  to  measure  simultaneously  the  real  and  imaginary 
parts  of  the  refi'active  index  of  a  spatially  sinusoidal  optical  grating  as  well  as  its  exact  location  in 
space  relative  to  the  beams  that  produced  it.  The  theory  and  experiments  are  outlined  in  paper 
3.9  below. 


Advanced  degrees  awarded  from  this  project  (12/1/93  to  10/31/94) 

5.1.  Ph.D.  degree  to  Ping  Xia.  See  paper  3.12  above. 

5.2.  Ph.D.  degree  to  Nansheng  Tang.  See  paper  3.13  above. 

Professional  personnel  supported  (12/1/93  to  10/31/94) 

6.1.  Dr.  R.W.  Hellwarth 

6.2.  Dr.  J.P.  Partanen 

6.3  .  Dr.  Ping  Xia 

6.4.  Dr.  Nansheng  Tang 

6.5.  Mr.  Huapeng  Guan 

Interactions  (12/1/93  to  10/31/94) 

The  interactions  by  the  members  of  this  project  are,  in  addition  to  the  talks, 
conferences,  workshops,  and  seminars  listed  in  Section  4,  the  consulting  services 
performed  by  Professor  Hellwarth  for  Lawrence  Livermore  National  Laboratory 
(contact:  Dr.  Mark  Henesian).  The  subject  matter  for  these  consultations  were  lasers 
and  laser  beam  interactions. 


MAJOR  ACCOMPLISHMENTS:  11/1/94  TO  10/31/95 


The  major  accomplishments  of  this  grant  period  were  as  follows. 

Ultra-high  electromagnetic  field 

During  the  current  period,  we  achieved  the  first  step  in  our  project  to 
determine  the  physical  limits  on  the  possible  strengths  of  electric  and  magnetic  fields. 
(The  highest  fields  imaginable  to  date  occur  in  extreme  astronomical  situations,  such 
as  a  neutron  star,  or  in  focused  laser  pulses.)  We  discovered  an  exact  focused-pulse 
solution  to  the  classical  MaxweU  equations  in  vacuum,  which  solution  gives  close  to,  if 
not  exactly,  the  strongest  electric  and  magnetic  fields  possible  for  a  ^ven  pulse  energy 
and  given  upper  limit  to  its  fi'equency  spectrum.  A  brief  account  of  this  achievement 
was  given  as  an  invited  paper  at  CLEO  '95.  (see  summary  reference  4.7  below  and 
Appendix  I.)  A  full  account  has  been  submitted  to  Phys.  Rev.  and  is  included  as 
Appendix  n.  We  are  currently  working  on  several  schemes  to  generate  these  solutions 
in  the  laboratory. 


Nonlinear  optical  properties  of  fullerenes. 

In  the  current  period,  we  achieved  more  measurements,  of  improved  accuracy, 
of  the  absorption  cross-section  (at  522  nm)  of  the  Cgo  molecule  in  its  lowest  excited 
singlet  and  triplet  states.  These  cross-sections  are  the  major  factors  governing  the 
excellent  performance  of  C^q  as  an  optical  limiter.  The  results  are  summarized  in  Pub. 
3.16  which  is  attached  as  Appendix  HI.  We  are  currently  extending  these 
measurements  to  determine  the  refractive  index  change  caused  by  excited  €50 
molecules. 


Atomic  vapors  as  nonlinear  elements. 

In  the  current  period,  we  achieved  the  first  direct  measurements  of  the 
diSusion  coefficients  of  excited  atoms  in  an  inert  buffer  gas.  We  measured  the 
diffijsion  of  both  cesium  and  potassium  atoms  in  their  lowest  excited  (4P)  levels  as 
they  dififiised  in  neon  and  xenon.  The  results  are  given  in  thesis  3.14  and  3.15  below 
and  in  manuscripts  being  prepared  for  publication.  A  brief  summary  is  given  in  the 
Digest  paper  4.9  which  is  attached  as  Appendix  TV. 

In  the  current  period  we  also  studied  a  trapped-atom  method  to  eliminate  the 
Doppler  broadening  which  tends  to  spoil  the  performance  of  atomic  vapors  as 
nonlinear  elements.  (See  Pubs.  3.4  and  3.6.)  Preliminary  results  were  presented  in  the 
Digest  Summary  (4.8  below)  which  is  attached  as  Appendix  V. 


Advanced  degrees  awarded  from  this  project  (11/1/94  to  10/31^5) 

5. 1  Ph.D.  degree  to  D.S.  GHassner.  See  paper  3. 14  above. 

5.2  Ph.D.  degree  to  Bing  Ai.  See  paper  3.15  above. 

Professional  personnel  supported  (11/1/94  to  10/31/95) 

6.1  Dr.  R.W.  HeUwarth 

6.2  Mr.  Huapeng  Guan 

6.3  Mr.  Bing  Ai 

6.4  Mr.  David  dassner 

Interactions  (11/1/94  to  10/31/95) 

The  interactions  by  the  members  of  this  project  are,  in  addition  to  the  talks, 
conferences,  workshops,  and  seminars  listed  in  Section  4.,  the  consulting  services 
performed  by  Professor  HeUwarth  for  Lawrence  Livermore  National  Laboratory 
(contact:  Dr.  Mark  Henesian).  The  subject  matter  for  these  consultations  were  lasers 
and  laser  beam  interactions. 


MAJOR  ACCOMPLISHMENTS:  11/1/95  TO  8/31/96 


The  major  accomplishments  of  this  grant  period  were  as  follows. 

Ultra-high  electromagnetic  field 

During  the  current  period,  we  continued  our  project  to  determine  the  physical 
limits  on  the  possible  strengths  of  electric  and  magnetic  fields.  (The  highest  fields 
imaginable  to  date  occur  in  extreme  astronomical  situations,  such  as  a  neutron  star,  or 
in  focused  laser  pulses.)  We  discovered  a  family  of  exact  focused-pulse  solutions  to 
the  classical  Maxwell  equations  in  vacuum,  which  solutions  give  among  the  strongest 
electric  and  magnetic  fields  yet  discovered  for  a  given  pulse  energy  and  given  upper 
limit  to  the  fi^equency  spectrum.  A  full  account  of  four  subclasses  in  this  family  was 
submitted  to  and  published  in  Physical  Review.  See  Publ.  3.18.  It  is  included  as 
Appendix  I.  We  are  currently  working  on  several  schemes  to  generate  these  solutions 
in  the  laboratory.  We  have  also  begun  a  collaboration  with  Professor  H.  Winfiil  of  the 
University  of  Michigan  Optical  Sciences  Laboratory  to  study  other  subclasses  of  these 
solutions. 


Nonlinear  optical  properties  of  fullerenes 

In  the  current  period,  we  achieved  more  measurements,  of  improved  accuracy, 
of  both  the  real  and  imaginary  parts  of  all  independent  elements  of  the  nonlinear 
susceptibility  that  governs  degenerate  four  wave-mixing  in  bulk  C50  over  the  entire 
wavelength  range  of  0.74  to  1.7|im.  We  observed  one  two-photon  resonance  and  its 
decay  to  a  clear  asymptotic  limit  at  long  wavelengths.  A  description  of  our  technique 
and  early  measurements  has  been  submitted  to  Physical  Review.  In  the  current  period, 
we  also  investigated  the  question  raised  by  the  well-established  five-fold  increase  in 
absorption  cross-section  of  C50  excited  to  its  lowest  triplet  state.  That  question  is 
whether  there  is  a  corresponding  change  in  the  optical  polarizability  that  accompanies 
the  increased  absorption.  We  find  that  indeed  the  excited  triplet  polarizability  is  —5% 
higher  than  the  ground  state.  Details  are  given  in  the  long  summary  of  talk  4. 1 1,  and 
included  as  Appendix  II. 

Optical  properties  of  protonic-doped  species 

We  have  continued  our  studies  of  the  nonlinear  optical  effects  of  various 
dopants  in  liquids  and  solids.  We  completed  a  study  of  protonic  doped  bis- 
benzothiazole,  a  chromophore  commonly  used  in  optical  polymers.  We  found  that  the 
absorption  cross-sections  of  the  first  excited  singlet  and  triplet  states  of  the 
chromophore  were  ~1.40%  and  1 1%  larger  than  that  of  the  ground  state  to  which  the 
chromophores  returned  in  ~  1.4  ns.  Details  of  these  results  are  found  in  Publ.  3  . 17  and 
attached  as  Appendix  III. 


New  four-wave  mixing  interferometer 

In  this  period  we  completed  development  of  an  interferometer  capable  of 
measuring  both  the  magnitude  and  phase  of  the  third-order  nonlinear  optical 
polarization  density  that  is  responsible  for  degenerate  four-wave  mixing.  We  used  this 
technique  successfuUy  to  find  the  complex  nonlinear  susceptibility  of  an  optical 
chromophore  dissolved  in  methylene  chloride.  Details  are  given  in  Publ.  3.19 


Experimental  demonstration  of  the  first  clear  example  of  a  strongly- 

coupled  large  polaron 

We  completed  experimental  and  theoretical  studies  of  the  mobility  of  electrons 
photoexcited  to  the  conduction  band  of  the  clear  insulating  n-type  cubic  sillenite 
crystal  Bii2SiO20  (n-BSO).  We  measured  the  mobility  over  the  temperature  range 
300  to  480  K.  This  was  sufiBcient  range  to  compare  with  the  detailed  predictions  well- 
established  of  polaron  theory  that  has  only  one  undetermined  parameter,  the  band 
effective  mass.  Taking  that  mass  to  be  2.0  ±  0.1  electron  masses  gives  essentially 
perfect  argument  with  our  results.  The  extensive  checks  to  establish  that  we  have  seen 
what  we  claim  above  are  complete  and  being  prepared  for  publication.  A  preliminary 
account  of  our  experimental  results  is  given  in  the  summary  published  with  Conference 
Talk  4. 10. 


3.  PAPERS  AND  Ph.D.  THESES  PUBLISHED  FROM  THIS  PROJECT 

(12/1/93  T04/31/97) 


3.1  "Bipolar  formation  and  nonlinear  optical  properties  in  bis-thienyl  polyenes,"  C.W. 
Spangler,  M.Q.  He,  J.  Laquindanum,  N.  Tang,  J.P.  Partanen,  and  R.W.  Hellwarth,  Mat.  Res. 
Soc.  Symp.  Proc.  Ser.,  vol.  138,  Electrical.  Optical  and  Magnetic  Properties  of  Organic  Solid 
State  Materials,  pp.  607-612,  December  1993. 

3.2  "Excited-state  absorption  cross-section  of  C50/'  N.  Tang,  H.  Guan,  J.P.  Partanen,  and 
R.W.  HeUwarth,  Proc.  S.P.I.E.,  vol.  2143,  pp.  272-280,  January  1994. 

3.3  "Studies  of  optical  nonlinearity  in  bis-thienyl  polyenes,"  N.Tang,  J.P.  Partanen,  and 
R.W.  HeUwarth,  J.  Laquindanum,  L.  Dalton,  M.Q.  He,  and  C.W.  Spangler,  S.P.I.E. 
Symposium  Proceedings  Series,  vol.  2285,  Nonlinear  Optical  Properties  of  Organic  Materials 
Vn.  pp.  186-195,  July  1994, 

3.4  "Optical  image  processing  by  an  atomic  vapour,"  I.  Biaggio,  J.P.  Partanen,  B.  Ai,  R.J. 
Kmze,  and  R.W.  Hellwarth,  Nature,  vol.  371,  pp.  318-320,  September  1994. 

3.5  "Nonlinear  optical  properties  of  C60  films  between  740  and  880  nm,"  F.P.  Strohkendl, 
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We  describe  and  analyze  an  all-optical  method  to  measure  both  the  magnitude  and  the  phase  of  each  (complex) 
element  of  the  tensor  that  determines  degenerate  four-wave  mixing  in  a  thin  isotropic  film  on  a  thick  sub¬ 
strate.  We  use  the  method  to  characterize  the  tensor  of  a  10-^tm  Ceo  film  (on  a  CaF2  substrate)  for  110-fs 
pulses  at  768  nm.  Using  a  fused-quartz  plate  as  the  nonlinear  standard,  we  find  the  two  independent  Maker- 
Terhune  elements  of  to  have  magnitudes  |ci22i(“^^j  “«*>)!  =  (0.44  ±  0.03)  x  10~^^  esu  and 

|cii22(-w,  0),  w,  -o))!  “  (0.50  ±  0.03)  X  10~^^  esu  and  phase  angles  <^1221  and  <^1122,  whose  magnitudes  are 
145°  ±  17°  and  139°  ±  10°.  ©  1997  Optical  Society  of  America.  [80740-3224(97)01901-2] 


1.  INTRODUCTION 

In  degenerate  four-wave  mixing  (DFWM),  three  beams  of 
the  same  frequency  co  overlap  in  a  medium  to  generate  ei¬ 
ther  a  fourth,  or  more,^  new  beams,  also  of  frequency  o). 
At  low  beam  intensities  the  power  in  each  newly  gener¬ 
ated  beam  is  proportional  to  the  product  of  the  powers  of 
the  three  incident  intensities,  and  the  constant  of  propor¬ 
tionality  has  as  a  factor  the  absolute  square  of  an  element 
of  a  third-order  susceptibility  tensor,  often  referred  to  as 
or  simply  As  the  beam  frequency  w  is  increased 
to  where  it  approaches  or  exceeds  either  (1)  the  lowest 
one-photon  resonance  frequency  of  the  medium  or  (2)  half 
of  the  lowest  two-photon  resonance  frequency  of  the  me¬ 
dium,  passes  from  being  a  positive  real  number  to  a 
complex  number  that  varies  in  a  complicated  way  with 
frequency  ct>.  Knowledge  of  the  phase  of  as  well  as  of 

its  magnitude,  is  relevant  for  the  understanding  of  the 
electronic  structure  of  materials  as  well  as  for  the  estima¬ 
tion  of  the  device  potential  of  a  given  material.  For  ex¬ 
ample,  some  devices  require  a  large  real  part  of  which 
results  in  field-induced  phase  changes,  whereas  other  de¬ 
vices  require  a  large  imaginary  part,  which  is  related  to 
field-induced  absorption. 

In  general,  the  third-order  susceptibility  is  defined 
as  a  function  of  four  frequency  arguments 
( —  0)4 ,  0)3 ,  0)2 ,  wi),  where  the  frequencies  o)i ,  o>2,  and  0)3 
are  the  frequencies  of  three  input  fields  that  create 
through  nonlinear  mixing  a  signal  field  with  frequency 
0)4  =  0)1  +  0)2  +  0)3.  For  example,  third-harmonic  gen¬ 
eration  is  related  to  3o),  o),  o),  o)),  whereas  DFWM 

is  determined  by  o),  o),  — o),  o)).  The  magnitudes 


and  the  phases  of  these  two  quantities  are  related  only  as 
o)  approaches  zero.  They  have  not  been  related  for  ex¬ 
periments  on  Ceo  to  date.  Phase  measurements  of 
with  third-harmonic  generation  based  on  the  interference 
of  substrate  and  thin-film  signals  have  been  frequently 
performed.  However,  most  DFWM  experiments  have 
measured  the  intensity  of  their  signal  beams  and  there¬ 
fore  only  the  magnitudes  of  various  tensor  elements. 
A  number  of  methods  that  measure  the  phase  and  the 
amplitude  of  ;^^^^(-a),  w,  -o),  w)  have  been  devised.  Ex¬ 
amples  are  the  DFWM  method  by  Tang,^  the  phase- 
conjugate  double  interferometer,^  measurement  of  the 
DFWM  signal  strength  of  a  substance  in  solution  versus 
its  concentration,^  Kerr-gate  heterodyne  detection,^  and 
the  z-scan  technique.®  For  a  given  thin  film  on  a  thick 
substrate,  none  of  these  methods  can  separate  the  contri¬ 
butions  from  the  substrate  and  the  thin  film  to  the  ob¬ 
served  nonlinear  signal.  Here  we  introduce  a  new 
method  that  allows  one  to  compare  both  the  magnitude 
and  the  phase  of  all  tensor  elements  of  a  given  thin 
film  to  those  of  a  thick  substrate.  As  an  application,  we 
characterize  fully  the  tensor  of  a  polycrystalline  Ceo 
film  on  a  cubic  CaF2  substrate  at  768  nm. 

The  basis  of  our  method  is  as  follows.  In  wave  mixing, 
the  wave  vector  of  electric  polarization-density 
wave  created  by  the  nonlinear  medium  response  may  or 
may  not  be  equal  to  the  wave  vector  of  an  optical  beam 
that  has  frequency  o)  and  satisfies  all  boundary  condi¬ 
tions.  If  it  is  not,  the  output  beam  generated  will  be  re¬ 
duced  in  power  by  a  phase-mismatch  factor  F 
=  [sin(x)/x]^  over  the  case  where  kj„ix  =  k^.  Here,  x 
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=  |kjnix  “  hJlJ2,  where  is  the  interaction  length  of 
the  beams  in  the  medium  a.  Our  method  compares  the 
power  in  a  generated  fully  phase-matched  beam  (F  ~  1) 
with  the  powers  of  beams  that  have  i^=l-^,  0<^ 
1 ,  for  the  thin  Cgo  film  but  F  ^  10“^  for  the  thick  CaF2 
substrate.  In  the  fully  phase-matched  beam  the  signal- 
field  contributions  of  film  and  substrate  interfere,  accord¬ 
ing  to  the  phase  difference  (f)  between  their  respective 
elements.  We  establish  the  interference  and  determine 
(f)  by  comparing  the  power  in  the  phase-matched  signal 
beam  (jP  =  1)  to  the  powers  in  signal  beams  that  origi¬ 
nate  (a)  solely  in  the  Ceo  film  {F  ~  1  -  S)  and  (b)  solely 
in  the  substrate  {F  =  1). 

In  the  following,  we  (1)  give  a  theoretical  discussion 
that  relates  the  observed  signal  fields  to  the  tensor 
components  of  thin  film  and  substrate  and  derives  their 
relative  phase,  (2)  demonstrate  our  method  experimen¬ 
tally  in  a  Cgo  film,  and  (3)  give  a  discussion  of  the  results 
that  leads,  because  of  redundancies  and  consistency 
checks  inherent  in  our  method,  to  corrections  of  the  long- 
pulse  x^^^  values  of  fused  silica  appropriate  for  short 
pulses  (—110  fs). 

2.  THEORY 

We  analyze  the  various  beams  radiated  at  frequency  w  by 
the  nonlinear  polarization  density  created  by  three  beams 
(of  frequency  co)  incident  upon  a  thin  film  on  a  thick  sub¬ 
strate.  The  z  axis  is  chosen  normal  to  the  sample.  The 
three  incident  beams,  called  1,  2,  and  3  in  Fig.  1,  are  ar¬ 
ranged  in  the  so-called  forward-four-wave  mixing  geom¬ 
etry  and  have  wave  vectors  (a,  a,  6^),  (-a,  a,  6^),  and 
(—a,  ~a,  6q,),  respectively.  Here  the  subscript  a  takes 
the  values  a,  f,  s ,  or  q  for  air,  film,  substrate,  or  fused 
quartz.  We  assume  that  a.  The  incident  beams 
give  rise  to  Re[P^^^(r)exp(-ZwZ)],  a  nonlinear  polarization 
density,  at  space-time  points  (r,  t)  that  generates  nine 
new  beams,  labeled  4  through  12  in  Fig.  1.  Beam  4  is  the 


Fig.  1.  Beam  geometry  for  degenerate  forward-four- wave  mix¬ 
ing.  Three  beams,  labeled  1,  2,  and  3,  copropagate  nearly  col- 
linearly  along  the  2  axis.  The  angle  of  intersection  in  air  be¬ 
tween  beams  1  and  2,  and  2  and  3,  is  denoted  by  6.  The  sample 
normal  is  parallel  to  the  z  axis.  The  three  incident  beams  inter¬ 
sect  inside  the  sample  and  generate,  because  of  third-order  non¬ 
linear  optical  interaction,  nine  signal  beams  that  are  labeled 

4- 12  in  a  plane,  which  we  call  the  observation  plane,  parallel  to 
the  sample.  Beam  4  is  the  phase-matched  signal  beam;  beams 

5- 12  are  non-phase-matched. 


only  signal  beam  that  is  perfectly  phase  matched.  By 
measuring  the  strengths  of  some  of  these  beams  (for  dif¬ 
ferent  combinations  of  incident  polarization  states)  rela¬ 
tive  to  a  beam  generated  by  the  substrate  alone,  and  rela¬ 
tive  to  a  standard  fused-silica  quartz  plate,  we  are  able  to 
deduce  both  the  magnitude  and  the  phase  of  each  element 
of  the  complex  fourth-rank  nonlinear  susceptibility  tensor 
-0),  a>,  (o,  -co)  of  an  isotropic  thin  film.  We  express 
the  total  optical  electric-field  vector  at  (r,  t)  in  the  film  as 
Re[E{,  exp(-iwZ)],  where  the  complex  electric-field  ampli¬ 
tude  is  understood  to  contain  any  spatial  dependencies 
such  as  the  superposition  of  several  beams,  i.e.,  E{^ 
^  E{,(r).  Then,  using  the  convention  of  Maker  and 
Terhune,*^  we  write 

CO,  CO,  -co)i  KKK  ,  (1) 

where  x^  is  a  complex  fourth-rank  tensor,  whose  elements 

ttf 

we  discuss  below.  When  the  frequency  arguments  of 
are  omitted,  they  are  assumed  to  be  (  — w,  co,  <o,  —co)  as 
in  Eq.  (1).  Note  that  in  the  Maker  and  Terhune  defini¬ 
tion  a  factor  (1/4)  on  the  right-hand  side  of  Eq.  (1)  that 
stems  from  the  multiplication  of  terms  of  the  form 
Re[E{,  exp(-fcoZ)]  is  absorbed  into  x^-  We  will  approxi¬ 
mate  the  total  complex  electric-field  amplitude  in  the  air, 
the  film,  the  substrate,  or  in  the  fused-silica  reference  by 

E:  =  S  E;  exp(ik;  ■  r),  (2) 

j 

where  the  sum  is  over  both  incident  and  generated  beams 
and  a  refers  to  the  medium  in  which  the  field  is  observed. 
In  the  absence  of  any  nonlinear  interaction,  Eq.  (2)  is  an 
exact  solution  of  Maxwell’s  equations  with  position- 
independent  Ej.  To  solve  Maxwell’s  equations  for  the 
generated  beams  4  to  12,  we  will  find  it  sufficient  to  use 
the  slowly  varying  envelope  approximation  or  SVEA  for 
Ef  (j  ==  4  to 7  =  12),  and  we  assume  that  the  incident- 
beam  amplitudes  EJ  {j  =  1  toy  =  3)  are  unaltered  by 
the  nonlinear  interaction. 

Consider  the  SVEA  equation  for  beam  5  (this  is  of  the 
most  general  form  we  will  encounter); 

dE^  12770)  ^  ^ 

IT  =  '  r“;E?E^E3“*  expii^klz),  (3) 

where  the  wave-vector  mismatch  A/e  5  follows  from  Ikg] 
=  n^co/c  and 

k^  =  kj  +  k^  -  k^  -  A^^z,  (4) 

which  is  required  by  boundary  conditions.  Here  is 
the  refractive  index  of  the  medium  a.  A  nonzero  wave- 
vector  mismatch  gives  rise  to  phase-mismatch  factors  F 
<  1  described  in  Section  1.  The  term  on  the  right-hand 
side  of  Eq.  (3)  is  the  only  term,  from  a  larger  number  of 
terms,  for  which  \\k^\  a.  The  other  terms  are  ne¬ 
glected. 

In  order  to  estimate  the  phase-mismatch  factor  F,  we 
derive  a  relation  between  the  external  beam  angle  in  air, 
0  (see  Fig.  1),  and  the  wave-vector  mismatch,  A^s.  The 
beam  angle  is 
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6  =  tan  ^  ,  "■=  2acfo).  (5) 

From  Eq.  (4)  one  sees  that  for  beam  5  the  wave-vector  in¬ 
side  a  medium  a  is  given  by  kg  =  (a,  3a,  65),  which  is 
determined  from 

=  lOa^  +  {blf  =  =  2a^  +  bl  = 

(6) 

We  approximate  6^,  by  ,  as  our  angle  of  incidence  is 
close  to  normal.  Comparing  Eq.  (6)  with  Eq.  (4),  one  sees 
that,  for  beam  5,  the  wave-vector  mismatch  is 

^kl  =  bl-  h^  =  4a^c/{n^o))  --  6^o)l(n^c).  (7) 

Note  also  that  in  our  symmetric  beam  geometry  the  wave- 
vector  mismatch  for  beams  5  and  6  is  equal. 

In  our  experiments  the  thickness  of  the  C^o  film  If  is  10 
yam,  the  refractive  index  nf=  2.0,^  0  =  2.86°,  and  \ 
~  768  nm,  so  that  the  phase-mismatch  angle  x 
=  A^{Zy2  =  0.05  and  the  phase-mismatch  factor  F  is 
well  approximated  by  1.  However,  if  we  recalculate  this 
angle  for  the  substrate  where  the  refractive  index  Ug 
=  1.43  and  the  thickness  Ig  is  1  mm,  we  find  A^|Zg/2 
7.13.  From  Eq.  (3)  we  see  immediately  that  the 
phase-mismatch  factor  F  reduces  the  substrate  signal  in 
beam  5  to  1%  of  what  its  value  would  be  were  A^g  =  0. 
We  varied  the  incident-beam  angle  6  with  the  sample  be¬ 
ing  the  pure  substrate  without  film  and  found  the  ratio  of 
the  power  generated  in  beam  5  to  the  power  generated  in 
beam  4  (for  which  the  wave-vector  mismatch,  A^,  is  iden¬ 
tically  zero)  to  be  predicted  to  within  experimental  accu¬ 
racy  using  Eq.  (7)  in  Eq.  (3), 

The  field  generated  in  phase-matched  beam  4  obeys 


dz 


12'7r<y  * 

- X^lE^E^Ef, 

n^c  1  j  2  » 


(8) 


where  a  =  f,  s,  or  q  for  the  film,  substrate,  or  quartz 
plate,  respectively. 

In  this  paper  we  will  study  only  the  signals  generated 
in  beam  4,  5,  or  6.  With  the  substrate  phase  mismatch 
calculated  above  and  with  the  values  that  we  obtain, 
we  find  it  consistent  to  neglect  the  substrate  contributions 
to  beams  5  and  6.  They  alter  these  beam  powers  by  less 
than  1%.  That  is,  we  assume  that  the  SVEA  equation  for 
beam  6  in  the  film  may  be  written  as 


dz 


12770; 

rifC 


(9) 


with  negligible  contribution  from  the  substrate. 

Note  the  differences  among  Eqs.  (3),  (8),  and  (9);  we  ex¬ 
ploit  these  to  advantage.  All  the  beam-overlap  correc¬ 
tions  are  the  same  for  Eqs.  (3),  (8),  and  (9),  making  com¬ 
parisons  among  their  corresponding  signal  strengths 
more  accurate  than  for  other  signal  beams  generated,  as 
beams  7  through  12  have  only  two  of  the  incident  beams 
in  the  right-hand  side  of  their  SVEA  equation.  Therefore 
beam-overlap  corrections  are  different  among  beams  7 
through  12,  and  beam  ratios  become  less  accurate  for  in¬ 
terpretation. 

When  all  beams  have  the  same  frequency  (o  in  isotropic 
media  (such  as  the  polycrystalline  €50  film  or  the  fused- 


quartz  reference  plate),  the  nonlinear  polarization  density 
(1)  has  the  particularly  simple  form 

=  3c?22iE:e„  .  E,  +  6c?i22E.E,  •  E*,  (10) 

where  again  a  =  s,  or  q  for  film,  substrate,  or  quartz 
plate.  The  complex  coefficients  C1221  ^1122  called 

the  two  independent  elements  of  the  tensor  in  Eqs.  (3), 

(8),  and  (9).  Their  frequency  arguments  are  understood 
to  be  (-W,  w,  a>,  — w).  Here,  following  the  notation  of 
Maker  and  Terhune,^  we  use  the  letter  c  instead  of  x 
the  tensor  components  of  X‘  When  all  incident-beam  po¬ 
larizations  are  parallel  so  that  =  xE,  a  polarization 
density  parallel  to  the  field  is  generated  by  Eq.  (10), 
whose  scalar  amplitude  is  often  written  3ciiiiE^E*,  i.e., 
^1111  ^1221  +  2c  1122-  Note  that  here  again  the  total 

field  amplitude,  E,  is  understood  to  contain  any  spatial 
dependencies  such  as  those  that  arise  from  the  superpo¬ 
sition  of  several  beams.  For  cubic  media  such  as  the 

CaF2  substrate,  x"^  has  three  independent  tensor  compo¬ 
nents.  We  will  not  elaborate  on  these  here  because  they 
have  never  been  measured  and  we  do  not  know  the  orien¬ 
tation  of  the  cube  axes  in  our  substrate.  The  CaF2  sub¬ 
strate  signals  will  serve  simply  as  an  intermediary  to  the 
normalization  of  the  Cgo  elements  to  those  of  fused 
quartz,  which  we  know. 

In  order  to  calculate  the  observed  signal  fields  in  air  af¬ 
ter  they  leave  the  sample,  we  correct  for  Fresnel  reflec¬ 
tions  by  assuming  that  the  electric-field  amplitude  of  any 
beam  is  multiplied  by  the  Fresnel  factor  in  passing 
from  a  medium  with  refractive  index  to  one  of  index 


t 


al3  - 


ria  + 


(11) 


Multiple  reflections  are  not  taken  into  account.  Negli¬ 
gible  errors  are  expected  from  this  approximation  since 
the  fraction  of  the  beam  power  reflected  from  the 
Cgo- substrate  interface  is  less  than  3%  and  since  the 
round-trip  time  in  the  10-yam-thick  film  is  longer  than  the 
pulse  duration.  Here  a,  will  take  the  values  a  (air),  f 
(film),  s  (substrate),  and  q  (fused  quartz). 

In  the  following  solutions  of  Eqs,  (8)-(10)  we  call  the 
beam  amplitudes  in  air  just  before  the  sample  Ej  (j 
=  1,  2,  3),  and  we  call  the  generated  beam  amplitudes  in 
air  after  the  sample  Ej^  (j  =  4,  5,  6).  Then  the  solution 
of  Eq,  (3)  gives  directly 


EqL  ~  ^  12'n'(wZy^[c^221®3  ®^2  ■  ^1122(^1^^2  ■ 

+  E2E1  ■  (12) 

The  solution  of  Eq.  (9)  gives 


EeL  “  Z127rft;Z^[c(22iEi  E2  •  E3  +  0(122(^2^3  *  ^i 

+  E3E2  •  (13) 

These  two  equations  relate  signals  in  beams  5  and  6  to 
the  magnitudes  of  the  independent  tensor  elements  of  the 
Ceo  film.  The  solution  of  Eq.  (8),  which  contains  film  and 
substrate  contributions,  gives 
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^4L  ~  ^  127rC0{[c{221®^2  ■*"  ^1122(^^1^3  ‘  ®2 

+  E3E1  • 

+  A^^iEiEsElZs^aZ-^^s^sa/^^s}-  (l'^) 

In  this  case  the  signal  power  in  beam  4  contains  an  inter¬ 
ference  term  between  the  complex  film  and  substrate 
terms.  We  obtain  the  magnitudes  of  the  substrate  co¬ 
efficients  by  observing  the  signal  in  beam  4  from  the  sub¬ 
strate  alone,  which  from  Eq.  (8)  has  amplitude 

E42^  =  i  1277 a>^ ®  :  EiE3E2^gi^g^sa/c Tig  ,  (15) 

The  ratio  of  this  signal  to  that  from  the  fused-quartz  plate 
alone  is  used  to  determine  the  substrate  From  Eq. 
(8)  the  fused-quartz  plate  alone  produces  a  field  in  beam 
4: 

®4L  ~  ^127rco[c^221®^2  ®1  *  ®3  ^1122(^1^3  ■  ^2 

+  E3E1  ■  E^)]lgtlgtga/cng .  (16) 

The  phase  angle  of  the  thin-film  tensor  component 
c{jki  relative  to  the  substrate  tensor  element  can  be 
derived  from  experiment  by  measuring  (a)  the  power  in 
beam  4  for  the  thin-film  sample,  S4(f,  s)  ^  \c{j^ 
+  yc^jkil^,  y  =  const.,  (b)  the  power  in  the  non-phase- 
matched  beam  p  for  the  thin-film  sample,  Spif ) 
a  \c{jki\^,  and  (c)  the  power  in  beam  4  generated  by  the 
substrate  only,  84(5)  oc  The  beam  number  p  can 

take  the  value  of  either  5  or  6  or  both,  depending  on  the 
polarization  state  of  beams  1,2,  and  3.  The  correct  as¬ 
signment  of  p  can  be  derived  from  Eqs.  (12)  and  (13)  by 
substituting  the  polarization  vectors  for  the  fields  and  re¬ 
quiring  that  the  same  thin-film  tensor  component  gener¬ 
ates  signals  in  beams  4  and  p.  In  the  case  of  and 
<^1122  ^  single  polarization  configuration  of  beams  1,  2, 


Table  1.  Powers  84,85,  and  85  Generated  in  Beams  4,  5,  and  6  of  Fig.  1  for  Various  Polarizations  of  the 

Input  Beams  1,  2,  and  3" 


Polarization 

Geometry 

a 

Ssif.s) 

SM 

b 

Sdf,s) 

SM 

c 

S4(S) 

SM 

d 

S4{f,s) 

a 

T 

1.34  ±  3% 

1.14  ±  5% 

7.22  X  10“2 

±  S% 

1.00  ±  ^ 

[% 

T 

1122 

1221 

xxyy 

1122' 

1122 

1122 

1122 

1122 

b 

t  t 

1.05  ±  4% 

1.46  ±  6% 

9.93  X  10-^ 

±  6% 

1.04  ±  6% 

1221 

1122' 

xxyy' 

1122' 

1122 

1122 

1122_ 

11221 

c 

T  T 

1.31  ±  2% 

1.29  ±  3% 

0.235  ±  - 

1% 

0.73  ±2% 

T 

1111 

Till' 

Till 

Till 

1111 

1111_ 

1111 

1111 

d 

-  T 

1.60  ±  7% 

1.65  ±  7% 

0.132  ±  i 

5% 

0.77  ±  8% 

T 

1122' 

1122 

xyyx 

1221 

1221 

1221 

1221 

1221 

“The  polarization  geometries  shown  give  the  polarizations  for  beams  1,  2,  and  3  in  the  same  order  as  these  beams  are  shown  in  the  observation  plane  in 
Fig.  1.  The  powers  are  normalized  to  the  power  S^iq)  of  beam  4  generated  by  the  fused-quartz  plate  for  the  same  polarizations.  The  arguments  (/*,  s)  and 
(s)  indicate  signals  from  film-on-substrate  and  substrate  alone.  The  rows  and  columns  of  the  table  are  each  labeled  by  a,  b,  c,  and  d  so  that  a  single  table 
entry  can  be  referred  to  elsewhere  by  matrix  notation;  e.g.,  (c,  b)  refers  to  the  entry  1.29  ±3%.  The  brackets  below  each  table  entry  contain  the  space 
indices  of  the  tensor  element  governing  the  generated  beam  written  over  the  space  indices  of  the  tensor  element  of  fused  quartz  responsible  for  the 
normalization  signal. 


and  3  suffices.  In  the  case  of  <^>1221  two  different  polariza¬ 
tion  configurations  are  required.  For  the  correct  assign¬ 
ment  ofp,  see  also  Table  1.  The  magnitude  of  the  phase 
angle  can  then  be  derived  from 

S,{f,s)  =  S^(f)  +  esds) 

+  2^^|Sp{f  )S4(s)  cos{<pijici),  (17) 

where  ^  ^  which  takes  in  our  case,  Ceo  on 

CaF2,  the  value  0.843. 

3.  EXPERIMENT 

Experiments  were  performed  on  a  polycrystalline  (10 
±  0.5)-/xm-thick  Cgo  film  on  a  1-mm-thick  cubic  CaF2 
substrate  of  unknown  orientation.  A  piece  of  fused  silica 
of  effective  thickness  (1.46  ±  0.07)  mm  served  as  a  refer¬ 
ence.  Beams  were  derived  from  a  mode-locked  femtosec¬ 
ond  Ti: sapphire  oscillator  with  a  highly  stable  amplifier 
operating  at  20  Hz,  which  has  been  described  elsewhere.^ 
The  center  wavelength  X  was  768  nm,  and  the  spectral 
width,  10  nm,  was  close  to  transform  limited.  The  pulse 
duration  (intensity  full  width  at  half-maximum,  FWHM) 
derived  from  second-harmonic  autocorrelations  at  the 
sample  location,  as  well  as  from  Fig.  2,  was  (110  ±  2)  fs, 
assuming  a  sech^  shape  for  the  temporal  intensity  profile 
of  the  pulses.  A  single  beam  was  split  into  three  beams 
of  similar  intensity.  These  beams,  beams  1,  2,  and  3  of 
Fig.  1,  passed  through  a  single  achromatic  lens  of  focal 
length  50  cm  such  that  they  intersected  in  the  focal  plane 
of  the  lens.  The  transverse  intensity  beam  profile  in  the 
focal  region  had  a  diameter  of  150  pm  (FWHM).  The 
beam  intersection  angle  d  is  2.9  ±  0,2°.  The  maximum 
total  intensity  at  the  sample  was  20  GW/cm^.  The  polar¬ 
ization  state  of  each  beam  was  controlled  with  zeroth- 
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order  quartz  half-wave  plates  and  subsequent  calcite  po¬ 
larizers.  The  polarizers  were  the  last  optical  elements 
passed  through  by  the  beams  before  hitting  the  sample  to 
ensure  high  polarization  purity  at  the  sample.  Beams  1 
and  3  could  be  independently  delayed  with  a  computer- 
controlled  delay  line.  The  temporal  beam  overlap  was 
optimized  with  two-beam  autocorrelation  signals  and  was 
within  ±5  fs  of  the  optimum.  Spatial  beam  overlap  was 
achieved  with  a  pinhole  of  50-yLtm  diameter.  An  indexed 
holder  for  the  sample  and  the  pinhole  established  sample 
coordinates  relative  to  the  pinhole,  which  allowed  (a)  put¬ 
ting  the  thin  film  in  the  plane  of  optimum  spatial  beam 
overlap  and  (b)  repeatable  sample  positioning.  All  ex¬ 
periments  were  performed  in  a  thin-film  region  of 
~0.3-mm  diameter.  An  aperture  in  the  observation 
plane  of  Fig.  1  is  used  to  select  single  signal  beams  for  de¬ 
tection.  This  aperture  is  large  enough  to  block  all  beams 
emerging  from  the  sample  and  can  be  translated  in  x  and 
y  directions.  A  lens  after  the  aperture  is  positioned  such 
that  it  images  the  interaction  region  in  the  sample  onto  a 
4  mm  X  4  mm  silicon  detector.  We  measure  all  beam 
powers  emerging  from  the  sample  simply  by  translating 
the  aperture  in  the  observation  plane. 

Often  when  one  mixes  three  simultaneous  write  pulses 
to  generate  a  fourth  signal  pulse,  the  delay  in  time  of  one 


Fig.  2.  Degenerate  four- wave  mixing  signals  for  all-parallel 
beam  polarizations  versus  beam  delay.  The  signals  in  beams  4, 
5,  and  6  emerging  from  the  thin-film  sample  (Ceo  film  on  a  CaF2 
substrate)  and  the  signals  in  beam  4  from  the  pure  CaF2  sub¬ 
strate  and  the  fused-silica  standard  are  shown.  Beam  3  was  de¬ 
layed  for  signals  in  beams  4  and  6.  Beam  2  was  delayed  for  sig¬ 
nals  in  beam  5.  All  data  sets  are  fit  to  Gaussians  to  guide  the 
eye.  Note  that  for  the  thin-film  sample  the  signal  in  phase- 
matched  beam  4  is  smaller  than  in  non-phase-matched  beams  5 
and  6.  This  demonstrates  destructive  interference  between 
thin-film  and  substrate  signals.  Note  also  that  the  thin-film  sig¬ 
nals  in  beams  5  and  6  are  equal,  as  expected  for  the  parallel  po¬ 
larization  configuration. 


of  the  write  pulses  still  allows  the  observation  of  the  sig¬ 
nal  pulse.  The  delayed  signal  arises  from  a  transient  re¬ 
mainder  of  an  initial  excitation  grating  produced  by  two 
earlier  write  pulses.  Much  has  been  learned  from  study¬ 
ing  these  transients.  However,  we  failed  to  detect  these 
transients  in  the  Cgo  film,  the  CaF2  substrate,  and  the 
fused-quartz  standard  at  our  working  intensity.  Typical 
signals  versus  the  time  of  delay  of  one  of  the  three  inci¬ 
dent  beams  are  shown  in  Fig.  2.  Note  that  the  Ceo  sig¬ 
nals  in  beams  5  and  6  are  equal,  as  expected  for  the  all¬ 
parallel  polarization  geometry.  The  Ceo  signal  in  strictly 
phase-matched  beam  4  is  smaller  than  in  non-phase- 
matched  beams  5  and  6.  This  clearly  indicates  destruc¬ 
tive  interference  between  the  thin-film  and  the  substrate 
signal  fields.  Note  that  the  absence  of  any  transient  sig¬ 
nals  makes  the  analysis  in  terms  of  a  tensor  appropri¬ 
ate. 

Each  measurement  series  (which  contained  measure¬ 
ments  of  signal  beams  4,  5,  and  6  for  the  thin-film  sample 
and  measurement  of  beam  4  for  a  film-free  region  of  the 
substrate)  started  and  ended  with  measurements  of  beam 
4  in  fused  silica.  The  fused-silica  measurements  were 
used  to  estimate  uncertainties  in  our  experiments  that 
were  due  to  instrument  drift  as  well  as  to  normalize  the 
values  of  the  thin  film  and  the  substrate.  To  obtain  a 
single  signal  trace,  we  averaged  each  data  point  shown  in 
Fig.  2  over  20  shots.  Acquisition  of  a  single  signal  trace 
takes  ~1  min.  To  ensure  that  we  were  measuring  the 
third-order  susceptibility,  we  checked  the  intensity  de¬ 
pendence  of  the  thin-film  signal  over  three  orders  of  mag¬ 
nitude  in  signal  strength  with  total  beam  intensities  at 
the  sample  up  to  twice  our  usual  working  intensity.  We 
confirmed  the  intensity  dependence  P  to  be  cubic  with  x 
=  3.0  ±  0.1.  At  five  times  our  working  intensity,  prob¬ 
ing  in  a  different  region  of  our  sample,  we  noticed  the  ap¬ 
pearance  of  delayed  signals  and  possibly  sample  damage 
detected  through  an  increase  in  the  signal  background. 
Care  was  taken  to  avoid  such  high-intensity  conditions  in 
our  experiments.  The  comparison  of  signals  from  the 
sample  region  in  which  we  performed  the  experiments  re¬ 
ported  here  to  signals  from  other  sample  regions  gave  no 
indication  of  any  sample  damage  that  was  due  to  beam 
exposure.  Our  results  for  thin-film  and  substrate  signal 
powers  normalized  by  the  corresponding  phase-matched 
fused-silica  reference  signal  are  shown  for  four  different 
polarization  geometries  in  Table  1. 


4.  EVALUATION  AND  DISCUSSION 

Using  the  16  data  numbers  in  Table  1,  we  calculate  here 
the  magnitudes  and  the  phases  for  all  nonzero  tensor 
elements  of  the  thin  polycrystalline  Cgo  film,  and  we  per¬ 
form  consistency  checks  made  possible  by  multiple  inde¬ 
pendent  measurements  of  C1122  ^1221  coefficients  of 

the  C0O  film.  In  this  process  we  will  use  the  fact  that  the 
c  1122  coefficient  of  the  thin  film  occurs  in  polarization  ge¬ 
ometries  a,  b,  and  d  of  Table  1  and  can  therefore  be  used 
to  determine  the  ratio  01122/^1221  ^^r  fused  silica  relevant 
for  short  pulses  (110  fs).  We  also  estimate  small  (3%) 
corrections  to  the  best  characterized  x^^^  tensor  element  of 
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fused  silica,  Ci221j  appropriate  for  our  short-pulse  condi¬ 
tions  under  which  some  Raman  vibrations  are  frozen  out. 
In  addition,  we  note  the  approximate  size  of  the  tensor 
elements  of  the  CaF2  substrate,  whose  orientation  is  un¬ 
known. 

Owyoung  measured  the  rotation  of  the  polarization  el¬ 
lipse  of  a  beam  of  694-nm,  10-ns  pulses  in  fused  quartz, 
relative  to  the  rotation  in  CS2.^^  This  rotation  is  propor¬ 
tional  to  Ci22i(  —  (o,  CO,  CO,  -co),  which  is  known  to  within  a 
percent  in  CS2  (through  accurate  determination  of  its 
Kerr  constant).  By  this  comparison  he  found 
Ci22i(~o)y  CO,  co,-co)  to  be  1.28  X  10“^^  esu  within  ±7% 
in  fused  quartz.  By  measuring  the  polarized  and  the  de¬ 
polarized  Raman  spectra  in  fused  quartz,  Hellwarth  et  al. 
determined  that  -^6%  of  this  value  arises  from  the 
nuclear  motions  in  fused  quartz. 

The  110-fs  pulses  we  have  employed  here  are  too  short 
to  be  affected  by  Raman  vibrations  below  200  cm“^.^^  A 
rough  estimate  of  this  freeze-out  effect  suggests  that  our 
pulses  feel  roughly  half  of  the  6%  Raman  contribution. 
(More  accuracy  is  not  warranted.)  Therefore  we  base  our 
estimates  of  the  absolute  values  of  c  1221  ^1122  poly- 
crystalline  Cqq  on  the  value 

CI221  =  1-24  X  10“^^  esu  (18) 

for  fused  quartz.  Corrections  for  freeze-out  effects  and 
the  ±7%  absolute  uncertainty  in  our  quartz  standard 
should  be  kept  in  mind  in  comparing  our  results  to  other 
experimental  results  on  Cgo-  The  entries  (a,  a),  (b,  b), 
(d,  a),  and  (d,  b)  of  Table  1  compare  the  two  independent 
components  of  x^^^  foi*  fused  quartz  with  the  C1122  coeffi¬ 
cient  of  C0O  •  From  these  we  determine 


c?i22/c?221  =  1-091  ±  0,013.  (19) 

This  is  13%  below  the  estimated  ratio  of  1.25  for  long 
pulses. For  the  purpose  of  estimating  the  absolute 
value  of Cjiii  for  Cgo?  we  combine  Eq.  (18)  with  Eq.  (19)  to 
obtain 

Clin  =  (3.95  ±  0.05)  x  10“^^  esu  (20) 

as  an  appropriate  reference  value  for  our  short  pulses. 

Using  the  standard  quartz  values  above,  we  obtain  the 
nonlinear  tensor  elements  for  polycrystalline  Cgo  that  ap¬ 
pear  in  Eq.  (10).  These  are  given  in  Table  2  along  with 
the  data  entries  of  Table  1  and  relevant  equations  that 
were  used  to  derive  them.  In  quoting  phase  angles,  we 
assume  that  the  tensor  elements  of  j^(“Ct),  co,  oj,  —00)  of 
both  fused  quartz  and  the  CaF2  substrate  are  real  and 
positive.  We  find  that  the  phase  angles  for  all  nonzero 
tensor  components  are  the  same  within  experimental  ac¬ 
curacy. 

Kajzar  et  al.  saw  two  pronounced  resonances  in  their 
measurements  of  third-harmonic  generation  versus  wave¬ 
length  (between  0.8  and  1.9  yu,m)  in  a  Ceo  film.  They  ob¬ 
served  a  broad  resonance  at  an  incident  wavelength  of  1.3 
jum,  which  they  suggested  is  a  two-photon  resonance  of 
the  character  that  would  produce  a  resonance  in 
—  0),  <w,  w,  — w),  also  at  1.3  /xm  with  a  width  of —0.2 
If  we  assume  that  this  resonance  dominates 
“CO,  o),  0),  —ft))  at  our  wavelength  of  768  nm,  then  a 
simple  Lorentzian  model  for  the  resonance  line  predicts  a 
phase  angle  for  our  x^^^  elements  near  160°,  declining  in 
value  if  a  positive  real  background  contribution  to  x^^^  is 


Table  2.  Results  for  the  Amplitudes  and  the  Phases  of  the  Third-Order  Nonlinear  Susceptibility  Tensor 
xH-o),  <0,  ft),  -ft))  of  Cgo  Obtained  by  Use  of  the  Data  of  Table  1  in  Eqs.  (12)-(16)  for  the  Optical  Fields 

Generated  in  Beams  4,  5,  and  6  in  Fig. 


Quantity 


Value 


kll22/Ci22ll 


I<^(l22l 


k^22ll 


l<^^22ll 


1.13  ±  0.05^ 

(1.42  ±  0.03)  X  10“^2 

(142°  ±  5°)^ 
(0.50  ±  0.01)  X  10"^2 
(0.49  ±  0.01)  X  10~^2  esu^ 

(138°  ±  10°)^ 
(140°  ±  12°)^ 
(0.44  ±  0.03)  X  10"i2  egyi 
(0.44  ±  0.01)  X  10“^2  eg^- 

(145°  ±  17°)* 


"Results  are  expressed  in  terms  of  the  two  independent  tensor  elements  0^221  and  0^122  defined  in  Eq.  (10).  The  parallel-polarization  element 
must  equal  0(221  ^^1122  •  ^*^ar  the  determination  of  the  phases  it  is  assumed  that  the  nonlinear  susceptibility  in  the  CaF2  substrate  has  a  zero 

imaginary  part.  In  the  other  footnotes  to  this  table  the  symbol  (a,  b)  refers  to  the  data  in  row  a  and  column  b  of  Table  1. 

^Average  of  the  ratios  of  (a,  a)  to  (a,  b)  and  of  (b,  b)  to  (b,  a),  weighted  by  the  inverse  squares  of  their  errors. 

‘^Average  of  data  entries  (c,  a)  and  (c,  b),  weighted  by  the  inverse  squares  of  their  errors  and  substituted  in  Eqs.  (12),  (13),  and  (16). 

‘^From  data  entry  (c,  d)  interpreted  by  Eq.  (14),  with  (c,  c)  giving  the  relative  magnitudes  of  cl^xx>  presumed  real  and  positive,  in  Eq.  (15). 

^Weighted  average  of  values  from  entries  (a,  a)  and  (b,  b). 

Walues  predicted  by  values  superscripted  (a)  and  (b)  above  and  assuming  <^(122  ~  4^1221  • 

^From  data  entries  (a,  d),  (a,  a),  and  (a,  c). 

^  From  data  entries  (b,  d),  (b,  b),  and  (b,  c). 

'Weighted  average  of  values  from  entries  (a,  b)  and  (b,  a). 

Value  predicted  by  values  superscripted  (a)  and  (b)  above. 

*  Value  obtained  from  entry  (d,  d)  in  which  1 0(221 1  is  obtained  from  the  average  of  entries  (d,  a)  and  (d,  b)  divided  by  the  square  of  the  ratio  superscripted 
(a)  above. 
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added.  Clearly  our  data  support  the  existence  of  this 
resonance. 

The  magnitude  of  the  Cmi  value  derived  by  us  for  the 
Cgo  film  is  the  smallest  value  reported  for  this  coefficient 
at  any  wavelength  and  a  factor  of  5  below  the  next  near¬ 
est  value  reported  by  Lindle  et  al.  at  1.06  /xm.^^  Our 
value  at  768  nm  being  lower  is  consistent  with  the  two- 
photon  resonance  hypothesis.  We  find  in  addition  that 
|cii22/ci22iI  =  1-13  ±  0.05,  which  we  have  not  yet  inter¬ 
preted. 

The  deviation  of  the  tensor  elements  of  the  cubic  CaF2 
substrate  from  isotropic  S3nnmetry  is  evident  in  column  c 
of  Table  1.  Since  we  do  not  know  the  orientation  of  the 
cube  axes,  we  cannot  reconstruct  this  tensor.  However, 
the  substrate  tensor  elements  appear  not  to  deviate  more 
than  —30%  from  the  values  of  polycrystalline  CaF2,  being 
—70%  of  those  of  fused  quartz.  (See  column  c  of  Table  1.) 


5.  CONCLUSION 

We  have  demonstrated  an  all-optical  method  that  allows 
the  measurement  of  the  amplitude  and  the  phase  of  all  in¬ 
dependent  elements  of  the  degenerate  four-wave-mixing 
tensor  of  an  isotropic  thin  film  on  a  thick  substrate.  The 
method  employs  the  interference  of  the  degenerate  four- 
wave-mixing  signal  from  the  substrate  with  that  from  the 
film,  normalized  to  other  signals  generated  by  the  same 
incident  beams.  Our  observation  of  a  negative  real  part 
in  every  element  of  at  768  nm  of  Ceo  supports  the  ex¬ 
istence  of  a  two-photon  resonance,  suggested  by  others  to 
exist  near  1.3  /xm. 
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-  (Received  10  May  1.996)  -  '  '  : 

We  determine  the  band  mobility  of  photoexcited  electrons  in  cubic  n-type  Bii2SiO20-  We  measure  a  . 
room-temperature  mobility  of  3.4  ±  0.5  cm^/(V  s)  that  decreases  monotonically  to  1.7  ±  0.3  cnr/(y  s)  . 
as  the  temperature  is  increased  to  200  We  show  that  electrons  in  Bii2SiO20  form  large  polarons.' 

Our  results  arc  predicted  by  strong  coupling  polaron  theory  if  the  band  mass  of  the  electrons  is  chosen 
to  be  2.0  ±  0.1  electron  masses.  We  determine  the  electron- phonon  coupling  constant  and  effective 
longitudinal  optical  phonon  frequency  required  for  this  prediction  from  the  available  infrared  reflectivity 
spectrum  of  Bii2SiO20-  [80031-9007(96)022 14-4] 


PACS  numbers:  63.20-Kr,  42.70 JMq,  72.10Di,  7Z20.Fr 

The  mobility  of  a  photoexcited  electron  drifting  in  the 
conduction  band  of  transparent  n-type  cubic  Bii2SiCho 
crystals  (n-BSO)  at  300  K  has  been  reported  to  lie 
between  3  and5  cm“/Vs  [1-3L  If  the  usualindependent- 
collision  model  were  used  to  describe  this  "drift,  the 
thermal  mean  free  path  (^0.2  nm)  would  be  less  than  the 
de  Broglie  wavelength  of  the  electron,  and  the  collision 
rate  times  the  Planck  constant  (the  uncertainty  in  the 
electron  energy)  would  be  an  order  of  magnitude  greater 
than  the  thermal  energy.  In  this  “strong  coupling”  case 
the  Boltzmann  equation  cannot  be  expected  to  apply.  We 
report  measurements  of  mobility  vs  temperature  in  n-BSO 
which  can  be  explained  well  by  a  “polaron”  theory  that 
avoids  the  limitations  of  the  Boltzmann  equation  [4-7]. 

An  electron  in  a  polar  crystal  polarizes  the  lattice  in  its 
neighborhood.  The  electron  moving  with  its  accompany¬ 
ing  lattice  distortion  is  called  a  polaron.  We  argue  that  ♦ 
a  photoelectron  in  n-BSO  constitutes  the  clearest  case  of 
a  strongly  coupled  “large”  polaron,  i.e.,  a  polaron  whose 
wave  function  extends  over  many  atoms,  so  that  the  elec¬ 
tron  can  be  thought  of  as  interacting  with  phonons  rather 
than  with  independent  atoms.  Previous  examples  of  large 
polarons  were  found  in  a  variety  of  alkali  halides,  where 
mobility  values  lie  above  12  cm*/V  s,  the  values  for  KCl 
and  KBr  at  330  °C  [8].  These  materials  are  closer  to  the 
case  of  independent  collisions.  Band  mobilities  lower  than 
in  n-BSO,  like  the  value  of  0.5  ±  O.i  cm“/V  s  found  in 
orthorhombic  KNbO;,  [9],  correspond  to  “small”  polarons. 
Small  polarons  have  a  wave  function  extending  less  than 
an  interatomic  distance  and  move  by  hopping  or  tunneling 
[10].  Many  even  smaller  mobility  values  are  reported  for 
various  insulators,  but  these  generally  reflect  the  effects  of 
shallow  traps  [U1 1]. 

n-BSO  has  a  sufficiently  large  linear  electro-optic 
effect,  so  that  charge  separations  can  be  easily  seen. 
Therefore,  we  determine  the  electron  mobility  using  the 
holographic  time  of  flight  (HTOF)  method  [9,11-14], 
in  which  two  interfering  laser  beams  excite  a  spatially 
sinusoidal  pattern  of  charge  carriers  in  the  bulk  of 
the  sample.  The  evolution  of  this  pattern  is  measured 
optically  by  observing  the  development  of  the  space- 


charge-induced  index  change  caused  by  the  separation 
of  the  free  carriers  from  their  excitation  point.  HTOF 
techniques  offer  two  distinct  advantages  over  regular 
time-of-flight  methods:  (1)  The  length  scale  is  set  by  the 
spatial  period  of  the  sinusoid  and  is  easily  varied.  (2)  The 
movement  of  the  photoelectrons  is  measured  optically, 
allowing  a  higher  time  resolution  (given  by  the  laser  pulse 
length,  which  is  30  ps  in  our  case).  HTOF  measurements 
can  be  performed  in  the  presence  of  an  applied  electric 
field  (drift  mode)  [13],  or  without  any  applied  field 
(diffusion  mode)  [9,14],  where  purely  thermal  diffusion 
of  the  photoexcited  charge  distribution  is  measured.  In 
this  work  we  use  the  HTOF  technique  in  diffusion  mode, 
which  has  a  number  of  additional  advantages:  (3)  no 
electrodes  are  needed,  (4)  no  particular  sample  shapes  are 
required,  and  (5)  no  uncertainties  are  introduced  because 
of  possible  internal  electric  field  variations  caused  by 
trapped  space  charge. 

The  experimental  four-wave-mixing  configuration  is 
shown  in  Fig.  1.  We  use  a  frequency-doubled  NdrYAG 
laser  that  produces  30  ps,  532  nm  pulses  at  a  repetition 
rate  of  5  Hz.  A  beam  splitter  sends  part  of  a  pulse 
into  a  delay  line,  to  act  as  a  probe  beam,  and  the 
other  part  is  split  again  into  two  write  pulses  that  arrive 
simultaneously  in  the  sample.  They  excite  a  sinusoidal 
electron  distribution,  ^sin(A"^,c),  from  donor  sites  with 
energy  levels  near  the  middle  of  the  3.2  eV  band  gap 
into  the  conduction  band.  Their  charge  is  compensated 
initially  by  that  of  the  newly  created  ionized  donor  sites. 
Diffusion  tends  to  make  the  conduction  band  electron 
distribution  uniform,  uncovering  the  space-charge  field 
Esc  of  photoionized  donor  sites  with  the  following  time 
dependence  [13,14]: 

Esciz,t)  ^  [1  Q\p(-t/T)]s\n(Kgz).  (1) 

where  =  tq  ‘  H-  tq  is  the  average  time  for 

photoexcited  electrons  to  remain  in  the  conduction  band 
before  going  to  uniformly  distributed  traps  of  unknown 
origin.  To  is  the  diffusion  time. 
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FIG.  1.  Diffracted  signal  vs  probe-pulse  delay  time  in  sample 
CT3,  at  a  grating  spacing  of  0.38  /x,m.  The  inset  shows  the 
experimental  arrangement.  Two  pulses  wl  and  w2  (532  nm, 
30  ps)  arrive  in  the  sample  at  the  same  time.  The  probe  pulse, 
which  is  counterpropagating  to  wl,  is  diffracted  by  the  space- 
charge  grating,  and  produces  a  signal  pulse  counterpropagating 
to  w2.  The  beam  splitter  BS  sends  the  signal  pulse  into  the 
detection  system.  The  beam  diameters  are  1 .7  mm  for  the 
write  beams  and  0.5  mm  for  the  probe  beam.  All  beams  are 
vertically  polarized.  The  grey  arrow  represents  a  cw  erasing 
beam  at  514  nm.  It  illuminates  the  crystal  from  the  top  with 
an  intensity  of  approximately  0.1  W/cm". 


Here  e  is  the  unit  charge,  is  the  magnitude  of  the 
sinusoidal  space-charge  modulation  wave  vector,  ytt  is 
the  band  mobility,  k  is  Boltzmann’s  constant,  and  T  is 
the  absolute  temperature.  This  analysis  assumes  that  a 
small  fraction  of  the  donors  are  photoexcited,  and  that 
the  evolution  of  the  photoelectron  pattern  is  dominated  by 
diffusion,  as  justified  experimentally  below. 

The  space-charge  field  of  (1)  modulates  the  refractive 
index  of  the  material  via  the  linear  electro-optic  (Pockels) 
effect.  We  detect  the  resulting  phase  grating  by  diffract¬ 
ing  the  time-delayed  probe  pulse  from  it  (see  Fig.  I).  We 
note  that  the  buildup  time  of  the  space-charge  field  (1) 
does  not  depend  on  the  write  fiuence  [14],  and  that  it  is 
given  by  the  diffusion  time  when  it  is  much  shorter  than 
the  free  carrier  lifetime.  We  know  that  photoexcitation 
of  holes  is  not  influencing  our  measurement,  for.  if  hole 
transport  were  significant,  we  would  not  observe  the  time 
dependence  in  ( 1 ). 

The  sample  is  homogeneously  illuminated  all  the  time 
by  a  514  nm  argon  ion  laser  beam  with  an  intensity  of 
approximately  0.1  W/cm“.  In  the  0.2  s  interval  between 
two  measurements,  this  illumination  erases  the  space- 
charge  induced  grating  created  by  the  write  pulses. 
During  the  measurement  time  (a  few  ns)  the  erase  beam 
deposits  6  orders  of  magnitude  less  energy  than  the  write 
pulses.  We  verified  that  it  does  not  affect  our  results. 

We  use  three  nominally  undoped  n-BSO  samples 
labeled  SUl,  CTl,  and  .CT3.  These  crystals  are  well 
characterized  by  many  experiments  as  described  in 


Refs.  [15,16].  The  SUl  sample  has  been  us^  in  previous- 
HTOF  experiments  in  the  drift  modc  |24  I43]rAK^ .  ^  • 
The  result  of  a  typical  mjeasurcmeht  ui"slk)w  &Hg.  1 
.  When  the  probe  pulse  precedes  the  write,  puls^^j^are' is 
onljra  small  signal.  A  relatively  strong  signal  caused  by"  : 
third  order  nonlinear  effects  is  observed  when  "the'  three 
pulses  are  present  in  the  crystal  at  the  same  thnel  When 
the  probe  pulse  is  delayed  (positive  times  in  the  figure), 
one  sees  the  diffracted  signal  increase  as  the  photoexcited  ' 
electrons  diffuse  away  from  the  positive  ions  at  which 
they  were  bound,  thus  creating  a  space-charge  field. 

The  signal  observed  before  zero  delay  has  amplitude 
Eq  and  is  due  to  scattered  light  from  the  probe  pulse  and 
to  the  remains  of  a  grating  that  is  not  completely  erased 
by  the  argon  laser  beam  (scattered  light  from  the  write 
pulses  was  taken  into  account  separately).  Although  it 
is  very  small,  this  background  grating  needs  to  be  taken 
into  account  when  measuring  the  buildup  time.  We  fit  the 
detected  diffraction  efficiency  77  to  7}(t)  ^  |Esc(0  + 
by  adjusting  £0  and  r  [see  (1)].  The  relative  phases  of 
Esc  and  £0  are  uncertain,  especially  at  the  low  intensities 
where  we  performed  the  buildup  time  measurements.  For 
each  data  set  we  perform  two  least  squares  fits:  once 
imposing  a  phase  shift  of  90'’  between  £0  and  and 
once  imposing  a  0  phase  shift.  To  obtain  a  final  value  for 
the  rise  time,  we  average  the  results  of  the  two  fits,  which 
differ  by  20%  at  most. 

In  order  to  determine  the  mobility  accurately,  we  must 
establish  that  we  are  in  the  low-fluence  limit  discussed 
above.  To  do  this  we  measure  the  fluence.  dependence 
of  both  the  buildup  time  and  the  magnitude  of  the 
signal.  Figure  2(a)  shows  the  grating  amplitude  and  the 
buildup  time  r  as  a  function  of  write  fluence  in  sample 
CT3.  For  small  fluences,  the  grating  amplitude  grows 
linearly  with  the  fluence  of  the  write  pulses,  and  the 
buildup  time  is  a  constant  independent  on  fluence.  As 
the  write  fluence  F  approaches  10  mJ/cm*  saturation 
effects  appear  [9].  From  the  data  in  Fig.  2(a)  one  derives 
a  photoexcitation  cross  section  a  between  and 

10“^^  cm*,  consistent  with  an  earlier  estimate  [15].  We 
performed  all  the  mobility  measurements  described  below 
with  a  write  fluence  of  the  order  of  I  mJ/cm". 

The  dependence  of  the  buildup  times  on  grating  spacing 
iTTjKa  is  shown  in  Fig.  2(b).  The  data  points  were  ob¬ 
tained  by  averaging  several  measurements.  The  buildup 
time  is  proportional  to  the  square  of  the  grating  spacing. 
This  confirms  that  the  electron  avoids  shallow  and  deep 
traps  for  so  long  that  the  buildup  time  r  corresponds  to 
the  electron  diffusion  time  to  and  the  grating  buildup 
is  dominated  by  diffusion  [see  (1),  r  =  tq].  The  elec¬ 
tron  lifetime  tq  in  the  SUl  sample  has  been  indepen¬ 
dently  determined  to  be  26  ±  2  ns  [13].  The  solid  line 
in  Fig.  2(b)  is  the  result  of  a  least-squares  fit  of  (2)  to  the 
data  of  both  samples,  with  only  the  mobility  as  a  tree  pa¬ 
rameter.  It  gives  a  mobility  value  of  3.4  ±  0.5  cm“/V s. 
The  buildup  times  in  the  SUl  and  CT3  samples  are  the 
same  within  the  experimental  error.  This  suggests  that  the 
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FIG.  2.  (a)  Grating  buildup  time  and  amplitude  vs  write  pulse 

fluence  in  sample  CT3  for  a  grating  spacing  of  0.56  The 
dashed  line  gives  a  linear  dependence  of  the  grating  amplitude 
on  fluence.  (b)  Log-log  plot  of  the  grating  buildup  time  r  vs 
grating  spacing  in  samples  SUl  and  CT3  at  1  mJ/cm“  fluence. 
The  solid  line  is  obtained  from  Eq.  (2)  using  a  mobility  value 
of  3,4  cm^/V  s. 


mobility  is  an  intrinsic  property  of  the  material.  This  idea 
is  also  supported  by  our  measurements  in  the  CTl  crys¬ 
tal  at  the  two  largest  grating  spacings,  where  we  found 
the  buildup  time  to  agree  with  the  results  in  Fig.  2(b). 
Our  mobility  value  is  consistent  with  other,  less  accu¬ 
rate  results  reported  previously  [1--3].  A  mobility  of 
50  cm'/V  s  derived  from  grating  buildup  time  in  a  similar 
experiment  [17]  is  erroneous  because  trapping,  not  diffu¬ 
sion,  dominated  the  dynamics  [18], 

Figure  3  shows  the  temperature  dependence  of  the  band 
mobility  in  CT3  and  SUl.  The  sample  was  enclosed  in  a 
temperature-controlled  oven  with  small  apertures  for  the 
laser  beams.  We  set  up  the  experiment  at  room  tem¬ 
perature  and  selected  an  appropriate  write  fluence  around 
I  mJ/cm“  as  described  above.  We  raised  the  temperature, 
stopping  every  20°-30°  to  stabilize  it  and  perform  a  mea¬ 
surement,  Care  was  taken  not  to  change  the  adjustment 
of  the  experimental  setup  when  changing  the  temperature 
of  the  crystal.  The  data  in  Fig.  3  were  obtained  at  a  grat¬ 
ing  spacing  of  0.38  ^am.  Other  measurements  performed 


FIG.  3.  Electron  band  mobility  vs  temperature  in  sample  SUl 
and  CT3,  and  predictions  of  Eq.  (4). 


at  longer  grating  spacings  from  0.43  to  0J7  fim  give  the 
same  temperature  dependence,  within  experimental  error. 

The  observed  decrease  of  the  mobility  with  rising 
temperature  and  the  polar  nature  of  the  n-BSO  lattice 
suggest  that  the  band  mobility  is  controlled  by  interaction 
with  longitudinal  optical  (LO)  phonons  [4-7].  The. 
interaction  of  an  electron  with  a  polarizable  lattice  has. 
been  studied  in  Refs,  [5-7]  starting  from  the  “Frohlich 
Hamiltonian”  which  has  three  material  parameters:  the 
electron-phonon  coupling  constant  a,  a  single  LO  phonon 
frequency  i'lo*  and  the  effective  mass  of  the  electron  in 
the  conduction  band.  In  the  alkali  halide  crystals  where 
much  of  the  polaron  theory  was  applied  [8],.  a  single  LO 
phonon  frequency  is  appropriate.  However,  /i-BSO  has 
many  polar  optical  vibrational  branches  with,  frequencies 
between  50  and  900  cm**  [19].  In  order  to  apply  the 
extensive  predictions  of  the  existing  polaron  models  [4- 
7]  to  our  results,  we  have  devised  two  mathematical 
schemes  to  imitate  the  phonon  structure  in  /i-BSO  by  a 
single  “average”  or  “effective”  LO  phonon  branch  having 
frequency  v\^q  Hz  and  oscillator  stren^  W  Hz.  These 
are  related  by  (W^/z^lo)  “  ^dc  -  Here  eao  =  5.7 

is  the  long  wavelength  limit  of  the  electronic  contribution 
to  the  refractive  index  squared,  and  =  50  is  the 
clamped  dc  dielectric  constant.  Our  schemes  determine 
I'LO  and  W  from  the  experimental  infrared  reflectivity 
spectrum  of  /z-BSO  at  room  temperature  [19].  In  terms 
of  J'LO  and  W,  the  dimensionless  coupling  constant  a 
defined  by  Frohlich  is 

or  =  (VK/z'lo)^  (3) 


where  Ry  =13.6  eV  is  one  Rydberg  of  energy,  h  is 
Planck’s  constant,  nC  is  the  effective  mass  of  an  electron 
in  the  conduction  band  (with  no  phonons),  and  rrie  is  the 
electron  mass.  Using  the  frequency  dependence  of  the 
dielectric  constant  of  n-BSO  we  obtain,  from  one  scheme, 
W/c  =  195  cm**  and  vlqIc  =  504  cm**,  and  thus  a 
coupling  constant  a  =  j rrie^  A  second  scheme 

gives  W/c  —  195  cm**  and  ploIc  =  495  cm**  instead, 
which  makes  negligible  difference  in  the  predictions. 

Since  our  experiment  was  performed  above  room 
temperature,  we  use  the  general  result  [Eqs.  (46)  and  (47)] 
of  Ref.  [7],  obtained  before  the  low  temperature  limit  was 
taken. 


e 


3v^  sinh(^/2) 


K{a,b). 


(4) 


where  ijl  is  the  mobility.  /3  =  hvi^olkT.  K(a,b)  = 
diiQO^{ii)/[Lr  +  or  -  bcos{vu)f^^,  =  (yS/2)^  * 
RyScoth(/3u/2),  b  -  R^/ s\nh{bl3v/2).  R  =  (v’  - 
w")/w'v,  and  v,w  are  temperature  dependent  varia¬ 
tional  parameters  [5.6].  We  find  that  putting  6  =  0  in  (4) 
makes  less  than  0.1%  error  throughout  our  temperature 
range.  This  is  useful  because  K{a,0)  =  K[{a)/a,  where 
K\  is  a  modified  Bessel  function  [20],  We  note  that  (4), 
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TABLE  1.  Variational  parameters  v  and  w  tor  various  cou¬ 
pling  constants  or  and  temperature  parameters  p. 


or 

V 

=  2 

w 

ar,= 

V 

3 

w 

a 

V 

4 

w 

1.5 

5.52 

429 

6.02 

4.01 

6.66 

3.70 

P  =  2.0 

4.48 

3.40 

4.94 

3.16- 

5.54 

2.88 

P  =  2.5 

3.89 

.  2.93 

4.32 

2.70 

4.90 

2.45 

where  K\{a)/a  can  be  substituted  for  K{a^b),  becomes 
equal  to  Eq.  (24)  in  Ref.  [21]  in  the  limit  of  small  a. 

We  determined  the  paraimeters  v  and  w  at  every  tem¬ 
perature  by  following  the  free  energy  minimization  proce¬ 
dure ‘described  in  Ref.  [6].  Table  I  gives  the  values  of  the 
v,  w  parameters  in  the  temperature  range  we  investigated. 
Using  these  values  we  estimate  a  polaron  radius  of  ap¬ 
proximately  0.6  nm  [22].  The  ;r-BSO  unit  cell  is  1.0  nm 
large  and  contains  66  atoms.  The  sphere  defined  by  the 
polaron  radius  contains  -^60  atoms.  Therefore,,  the  con¬ 
tinuum  approximation  of  Refs.  [5-7]  can  be  applied. 

The  only  unknown  parameter  in  (3)  and  (4)  is  the 
electron  effective  band  mass  m*.  The  prediction  of  (4) 
at  r  =  300  K  corresponds  to  our  room-temperature  mo¬ 
bility  value  of  3.4  cm-/Vs  when  setting  m*  =  2.01m^,. 
From  m*  and  (4)  we  can  predict  the  temperature  depen¬ 
dence  of  the  mobility.  The  result  is  shown  in  Fig.  3  to¬ 
gether  with  two  other  curves  obtained  using  =  \,lme 
and  m*  =  23me,  The  agreement  with  the  m*  = 
curve  is  very  good.  No  parameter  was  adjusted  to  fit 
the  experimental  temperature  dependence,  which  is  given 
mostly  by  the  effective  phonon  frequency  p\^q. 

In  conclusion,  we  have  presented  comparative  and  ther¬ 
mal  evidence  that  photoelectrons  in  nominally  undoped, 
n-type  Bii^SiO^o  are  large,  strongly  coupled  polarons; 
their  observed  band  mobility  is  intrinsic.  Their  small  mo¬ 
bility  makes  photoelectrons  in  n-type  Bii2Si02o  the  clear¬ 
est  example  of  such  polarons. 
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